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Abstract 

We study slip boundary conditions for simple fluids at surfaces with nanoscale chemical hetero- 
geneities. Using a perturbative approach, we examine the flow of a Newtonian fluid far from a 
surface described by a heterogeneous Navier slip boundary condition. In the far-field, we obtain 
expressions for an effective slip boundary condition in certain limiting cases. These expressions 
are compared to numerical solutions which show they work well when applied in the appropriate 
limits. The implications for experimental measurements and for the design of surfaces that exhibit 
large slip lengths are discussed. 
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I. INTRODUCTION 



The no-slip boundary condition was considered to have been experimentally established 
for simple liquids in the early 20th century. However, the refinement of a number of mea- 
surement techniques has recently led to the observation of nanoscale, and even micron-scale, 
violatic)ns of the no-slip boundary condition by simple fluids flowing over non-wetting sur- 
faces pj]. In many instances however, poorly controlled microscopic factors that influence 
the measured macroscopic slip length, such as roughness, chemical heterogeneity, or con- 
taminants such as air bubbles, have lead to apparent discrepancies in the magnitude of slip 
reported in the literature ^, 3, 4]. Thus, it is important to distinguish between effective or 
apparent slip, typically measured in macroscopic experiments, which emerges from the in- 
teraction of microscopic chemical heterogeneity, roughness and contaminants, and intrinsic 
slip, which arises solely from the chemical interaction between the liquid and a homogeneous, 
atomically flat surface. 

Slip is usually described in fluid mechanics by the Navier slip boundary condition Q]. 
This states that at a solid boundary, z = 0, the slip velocity, u, is proportional to the shear 
rate, dzU i.e. 

S d,u\^^o = (1) 

where the constant of proprtionality S is called the slip length. In some instances, experi- 
ments have found that the slip length can range from nanometers 6!l to tens of micrometers 

Q 

[7|. As slip on this scale can profoundly affect flows in micro and nanofluidic devices, these 
findings have generated considerable interest 8|]. For instance, large effective slip lengths 
potentially offer new ways of controlling flows in microdevices jol, [lo|. From a theoretical 
point of view, neither intrinsic nor effective slip lengths can yet be predicted microscopically. 
Nonetheless, a useful way to study intrinsic slip is through atomistic computer simulation, 
using techniques such as molecular dynamics. Such studies suggest that flows over flat hy- 
drophilic surfaces will exhibit intrinsic slip lengths less than a few nanometers^jwhile flows 
over flat hydrophobic surfaces should have slip lengths of tens of nanometers pjj. Indeed, 
strong experimental support for this picture is now beginning to emerge a]. 

However, it is on so-called superhydrophobic surfaces that slip lengths as large as tens of 



micrometers have been observed 



(see flgured]). The best known example of a superhy- 



drophobic surface is the leaf of the lotus plant, which possesses a microstructure and surface 
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FIG. 1: The figure illustrates the flow over a highly heterogeneous hydrophobic surface charac- 
terized by length scale L. At low pressures, the liquid is in the Cassie state, where it does not 
penetrate into the surface, leading to large effective slip lengths (left). At intermediate pressures 
the liquid begins to penetrate the surface (described by radius of curvature R) and the effective 
slip decreases (center). Finally, at sufficiently high pressures the liquid penetrates the surfaces, 
which will drastically reduce the effective slip length (right). 



chemistry that prevents water from wetting its surface, leading to droplet contact angles 
close to 180" 13|. Recently, nanotechnologists have learnt to mimic this so-called Lotus 

n 

effect by creating superhydrophobic surfaces [M| using materials such as carbon nanotubes 



assembled in dense forests 151]. The repulsion of water by such surfaces means that droplets 
or larger scale flows are essentially lubricated by a layer of air, leading to what is clearly a 
large effective slip length, with drag only occurring at the few points of the surface where 
the flow makes contact with the substrate. Again there is no rigorous theoretical description 
of how such effective slip lengths depend on the underlying microstructure of such highly 
heterogeneous surfaces. 

It is of interest then to study how effective slip lengths emerge from heterogeneous intrin- 
sic slip lengths. Such problems have been studied by numerical methods including molecular 
dynamics [3], lattice Boltzmann simulations [l7] and numerical solutions of the Stokes equa- 
tions 13]. In addition, some exact solutions are known for flow in channels both over and 



along stripes of alternating no-slip (5 = 0) and perfect slip {6 = oo) [19|, l20|, l2l|. However 
we are still unable to deduce effective slip lengths from the underlying microstructure and 
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chemistry of a surface in the general case. A deeper understanding of effective shp would 
give insight into how artifacts such as roughness and nanobubbles effect experimental mea- 
surements of slip. It may also allow optimization of superhydrophobic surfaces to extremize 
slip lengths for use in devices jo], llo| . 

In this paper our goal is to calculate effective, far-field slip lengths on chemically hetero- 
geneous surfaces, including nanoporous surfaces and surfaces covered in nanobubbles, which 
possess finite slip lengths < 5 < oo. In contrast, previous work has generally focused on 
the case where 5 is or oo. We begin by considering the hierarchy of length scales present 
in the problem and use this to define several distinct sets of problems that arise from this 
hierarchy in limiting cases. We then derive approximate solutions to two these problems 
to arrive at expressions for the effective slip length in the corresponding limiting cases. Fi- 
nally, we discuss the implications of these expressions for experiments and for the design of 
surfaces with large effective slip lengths. 



II. ANALYSIS 



The general problem we will examine here concerns the Stokes flow of an incompressible 
fluid past a surface described by a flnite slip length that is a function of position on the 
surface. We will consider simple shear flows, as might be encountered in many experimental 
situations. Thus at some distance W away from the heterogeneous surfaces, we apply a 
shear in the x-direction either at constant rate 7 or with a constant velocity u = Ug in the 
x-direction. For the constant shear rate, this results in the following boundary conditions 
at z = W : d zu\^^y^ — 7 where u is the x-component of the velocity u. For a shear flow 
past a flat surface which lies in the xy plane, z = 0, the slip length 5 = S{x,y) leads to the 
boundary condition S{x,y) dzu\^^^ = u\z=o. Note that we will initially ignore the effects of 
roughness induced by curvature of the liquid- vapor interface jisl as shown in flgure [H This 
amounts to assuming that the radius of curvature of the interface R is much larger than the 
other length scales in the problem. The effective slip length is then given in the limit as 

— > 00 by the expression 

6eff 5^mL=oo = m|^=oo. (2) 

There are three cases we have considered: 6 = 6{y) where 5{y) is periodic with period L 
(i.e. the shear is parallel to the patterning), 6 = 6{x) where 6{x) is periodic with period L 
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(i.e. the shear is perpendicular to the patterning) and S = S{x, y) where S{x, y) is periodic in 
the X and y direction with unit cell [L^, Ly) (i.e. the flow occurs over rectangular patches). 
We will focus on slip lengths which are patterned in stripes or patches with sharp edges so 
that S{x,y) will generally be considered to be a piecewise constant function. We will also 
assume that < S{x,y) < oo. For instance, in the case where the stripes oriented parallel 
to the direction of shear: 

^ Si 0<y<a 
62 a < y < L 

where a is the stripe width and Si < 62- In figure [H 61 would be the slip length of the solid 
surface, 6s and 62 would be the slip length over the vapor regions, 6g. In what follows we 
will restrict ourselves to presenting the analysis for the first geometry, where S = S{y), and 
simply report the analogous results of our calculations for the other two cases. Although 
the analysis is simplest for this first case, the approach in latter two geometries does not 
differ significantly from that presented here. 

It is useful at this stage to consider the magnitudes of the relevant length scales in the 
problem. Here we will assume that intrinsic slip lengths for smooth solid surfaces are at most 
10-20 nm, consistent both with recent measurements of the the slip lengths for hydrophobic 



surfaces jsl and with the results of molecular dynamics simulations 



b estimate slip 



lengths at the liquid- vapor interface, we will use De Gennes's expression 22|: 5^ ~ {fii/fig)t 
where is the viscosity of the liquid (gas) and t is the thickness of the gas layer. For pure 
water flowing over air at room temperature, we estimate that 6g ~ 50t. The length scales 
that describe the patterning can clearly vary widely. For the thiol functionalized carbon 
nanotube forests studied in Ref 12|], L ~ 1 fim, a ~ 50 — 100 nm and t > 1 fim, giving the 
following hierarchy of lengths: 6s <^ a <^ L <^ 6g. The mixing device _*ed in Q 
has grooves of width L — a ~ 100 fim and depth t ~ 50 fim, spaced at L ~ 10 mm, giving 
6s {L — a) <^ 6g ^ L. In contrast, a hydrophobic surface contaminated by nanometer- 
sized bubbles or a hydrophobic nanoporous surface might be described by the hierarchy: 
a L <^ 6s and L <^ 6n. The nanostructured hydrophobic channels studied by molecular 

nn 

dynamics simulations in Refs. pj|, 118| would also be likely to satisfy this hierarchy. More 
typically however, unless it has been specially prepared or contaminated by vapor bubbles, 
we might expect a heterogeneous surface to be described by (5i ~ ^2 <^ a ~ 

The hierarchy of length scales present will effect the magnitude of the terms in the slip 
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boundary condition. Introducing non-dimensional spatial coordinates x = x/L, y = y/L 
and z = z/L, the boundary condition for stripes oriented parallel to the direction of shear: 



(Si/L) diu\^^Q = u\b=o, <y < P 

{62/L) d^u\^^^ = M|i=o, [3 <y <1, (3) 

where j3 = a/L. In what we expect to be the most common situation, where 5i and 82 <^ L, 
the boundary conditions are no-slip to zeroth order: 

n|,=o = 0(5i/L), 0<y</3 

M|i=o = 0((52/L), I3<y<l. (4) 

For hydrophobic surface covered by nanometer scale bubbles, or a nanoporous substrate 
with L < 10 nm, we would have L -C 5^ and L <^ 6g, so that the boundary conditions are 
shear free at zeroth order: 

dM,=o = 0{L/6s), 0<y<(3 

d,u\,=o = 0{L/6g), P<y<l. (5) 



For superhydrophobic surfaces, such as those fabricated in Ref . 12] , which satisfy 5^ ^ L ^ 
6g, the boundary conditions are mixed at zeroth order: 

uU=o = Oi6jL), 0<y<P 
d,u\,^, = OiL/Sg), P<y<l. (6) 

Surfaces can be similarly defined both in the case of stripes perpendicular to the shear 
6 = S{x), and in the case of patches or more complex patterns where 6 = 6{x,y). In these 
cases we consider (3 to be the area fraction of the solid (or more generally the area fraction 
of the surface with the smaller slip length ^i). 

Our approach is to treat these problems perturbatively, solving them to first order in 
the relevant small parameters. As we will show below, this perturbative approach succeeds 
for boundary conditions (jlj) and ([5]) but fails for boundary condition ([6]). In addition, 
for surfaces that satisfy (E]), exact solutions of the Stokes equations with the zeroth order 
boundary condition are only known in the case where the shear is parallel or perpendicular 
to the stripes 3]. 
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In fact, when the shear is parallel to the stripes {S = S{y)), as both the velocity and the 
pressure are a function of y and z only, the equation for the x-component of the velocity, u, 
is just Laplace's equation: 

diu + dju = 0. (7) 

Indeed this is the easiest geometry to treat, and as such we will generally restrict ourselves 
to presenting the analysis of this particular case. Nonetheless, apart from case we have 
been able to extend our calculations of effective slip lengths to the more general patterning 
where 6 = 6{x, y). 

We now seek a solution of ([7]) that satisifies the boundary condition ([5]) of the form: 



u 



[x,z) = uo{y,z) + aui{y,z) + O (a"^^ (8) 



i.e. an asymptotic series in a = L/5i. As equation ([7]) is linear, each of the terms in the 
expansion, Ui, are solutions of ([7]). The z = boundary condition ([5]) at zeroth order in a 
is just a shear-free condition: 

^(y.o) = o (9) 

for < y < 1 and the solution at this order is simply a homogeneous shear-free flow with 
Mo = Us- At first order in a the slip boundary condition becomes 

dBUi{x,0) = Us <y < P 
d,uiix,0) = P <y <1. (10) 

Further conditions follow from the lack of a pressure head, the periodicity of the flow in the 
y-direction (period L) and the fact that the velocity component of the flow normal to the 
surfaces a.t z = and z = oo must vanish. 

To solve the first-order problem, we use the periodicity of the flow in the ^-direction to 
write Ui as a Fourier series as follows: 

oo 

ui{y,z) = '^Un{z)exp{ikny) (11) 

n=0 

where kn = 27m and 

Un{z) = [ ui{y,z)exp{ikny)dy. (12) 
Jo 

Inserting ( ITTl) into (jZl) we find that 

Uo {z) = Ao + Boz (13) 



and 

(z) = A^e-'^"^ + 5„e'"^ (14) 

for n > 0. Furthermore, for n > as the upper boundary conditions on Un (z) apply at 
z = W/L ^ 1, the coefficients i?„ will be of order exp {—knW/L), so we may neglect these 
in the far field. 

The slip condition at z = gives 

dUo dui 
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% 



u. ( P + [1 - H)^^] (15) 



SO that 



f/o (z) =us{z- W/L) ([3+ {I- I3)^£\ . (16) 



Thus the solution to the first-order problem is given by 

It is easily verified then that to first order in 5s/ L, as u (x, z) = Ug + aui + O (a^), the 
effective slip length for fiow over parallel stripes is given by 

' (18) 



Seff 5i 62 

The analysis in the remaining two geometries is similar {6 = 6{x) and 6 = 6{x,y)), 
although somewhat more complicated as the Stokes equations do not reduce to the Laplace 
equation (equation ([7j)) in these cases. Nonetheless, both in the case of shear directed 
perpendicular to stripes of fractional width /? and in the case of regular patches of area 
fraction (3, we again find that this relation holds. We note that these relationships have 



previously been o 
Stokes equations 



Dserved to hold empirically for numerical solutions of the steady state 



isl . In the second case, ([5j), where and 62 are much less than L, a 
similar analysis to that given above reveals that 6eff = (35i + (1 — 13)52- Furthermore, if 
^1 ~ ^2 ^ L, equation (ITS]) also reduces to 5eff — [35i + {1 — (3) 82- 

These expressions can be tested numerically. Figure [2] compares the effective slip length 
inferred from by finite difference solutions of the Laplace equation for the case of flow directed 
along stripes {5 = 8{y)) with equation (fTS!) as a function of stripe width P for 5i = 0.1 — lOL 
and 62 = lO^L. The plot shows that equation (fT8l) gives an excellent approximation to the 
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FIG. 2: Effective slip length as a function of stripe width f3 as given by equation (jlSp (solid line) 
and by finite difference solutions of the Laplace equation (symbols) for the case of flow over parallel 
stripes for = 0.1 - lOL and 82 = lO^L. 

effective slip length for 61 > l.OL but underestimates the effective slip length by up to a 
factor of 2 for (5i = O.IL while still tracking the dependence of the slip length on (3. Similarly, 
Figure [3] shows the effective slip length as a function of 62 for 61 = 0.1 — lOL with (3 = 0.1 for 
flow parallel to stripes. Again it is apparent that equation (ITSl) is very successful for 61 > L 
but underestimates 6eff for 61 = O.IL as one might expect. 

Unfortunately, as noted above, the perturbative approach used here fails for the third 
case (EI), where 5s ^ L ^ 5g. The exact solution to the zeroth order problem is only 
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FIG. 3: Effective slip lengtli for f3 = 0.1 as a function of 62 for several values of 5i as given by 
equation (fTH]) (sold line) and by finite difference solutions of the Laplace equation (symbols) for 
the case of flow over parallel stripes. 

known for the parallel stripes case {6 = S{y)), and this solution is not differentiable at the 
heterogeneous surface [lol causing the expansion ([8]) to fail. At zeroth order (in 6s/ L and 
L/Sg), for the parallel stripes, the effective slip length is given by 6^jj- = ^ In (esc /57r/2) 



19(1 . Note that there are also solutions known to this problem in the transverse stripe case 



{6 = S{x)) [19|], which differs by a factor of one half from the perpendicular case, but not to 
the more general case (5 = 6{x, y)). 

Finally, in this section, we note that ( fT8l) can be generalized to more complex patternings 
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S = S{x,y). Provided S{x,y) ^ L everywhere in the domain then a similar analysis to that 
given above yields 

where the angle-brackets denote that mean value of the function over the surface. Similarly, 
if S{x,y) <^ L everywhere on the surface 6eff = {6{x,y)). Thus the results presented here 
are not restricted to sharp, well-defined patternings. 

III. DISCUSSION 

The implications of (ITS!) for the effective slip lengths of hydrophobic surfaces contaminated 
by nanobubbles induced by roughness such as the superhydrophobic surfaces studied in 
Ref. pJl] are interesting. Here, if a ~ t ~ L < 10 nm for instance, we would expect 6g ~ 
50 if: ^ 5s ~ 10-20 nm > L, so that according to our results here, 6eff ~ Ss/p. Thus 
the apparent increase in slip length for a hydrophobic surface contaminated by nanoscale 
bubbles remains proportional to 6s and will be effectively independent of 6g. A surface with 
effective slip described by f|T8l) could exhibit slip lengths several times larger than 6s e.g. for 
a/L = 0.5, /5 = 0.25 so that 6eff could be as large as 40-80 nm. This is certainly consistent 



with many experimental measurements of slip (e.g. see the review Ref. {4]). 

It is important to note however that we have not considered the effects of roughness 
in our calculations above. Following Ref. |23|, one can consider the Navier slip boundary 
condition for fluid flow past a general solid surface: 

t ■ ^\z=z(x,y) = ^ti^j ( ViMj + '^jUi)\,^,(^,^y) (20) 

where t and n are the tangent and normal to the surface respectively. This can be rewritten 
in terms of the radius of curvature of the surface, R, at each point on the surface as follows: 

t ■ U|^=^(^,j^) = _ ^^^^ n ■ V (t ■ U) \z=z{x,y) (21) 

as originally noted in Ref. |23|. Thus, by using (12T|) to incorporate the effects of curvature 
of the liquid-vapor interface on the slip length at this surface, the effective slip length can 
be written 
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From this expression it can be seen that the curvature will become important when i? ~ 5^ 
and will decrease the effective slip length. 

As noted above, when Sg < L, equation fllSp fails as can be seen in the numerical cal- 
culation (figures [2] and [3]). For parallel or perpendicular stripes, when 5^ -C L ^ 6g, 
we expect that 5e// should be given by the expressions due to Philip [l^: that is 
6eff ^ L In (esc (3tt/2). Although neither exact nor approximate analytic solutions are known 
for the two-dimensional case, experiments and numerical solutions [2^ suggest that 
6eff ^ L also for fixed f3 in this limit. Thus for fixed (3 there appears to be a cross-over from 
6eff ~ Ss to 6eff ^ L as 6s goes from above to below L. We note that Ybert et al. [241 have 
suggested that the expression: 

^eff ~ (23) 

may approximately interpolate between these two limits. This is consistent with ([6]), which 
would lead one to expect a correction to the zeroth order solution due to the finite slip length 
of the solid proportional to 6s. 

Finally, these results suggest that very large slip lengths (> lOO's of nanometers) cannot 
be achieved by structuring a hydrophobic surface on length scales of lO's of nanometers as 
6eff ~ 6s in this case. However, it suggests that a hierarchy of length scales, which can lead 
to considerable enhancements in contact angle, could also be a useful way of maximizing 
effective slip length. If a hydrophobic substrate were patterned both on nanometer length 
scales L' < 6s and on micrometer length scales L, one might enhance 6'^jj^ for the solid by a 
factor of 3-4 as discussed above. According to equation (1231) . this could lead to a substantial 
increase in the overall effective slip length 6eff if 6'^jj- ~ a or larger. 



IV. CONCLUSION 



In summary, we have considered surfaces with alternating stripes or patches of slip length 
61 and 62 patterned on a length scale L. In the far-field, we derived expressions for the 
effective or apparent slip length in several cases. When 61 and 62 L or when 52^*^1^ L, 
the effective slip length is the area weighted average of the two slip lengths: 6eff = P6i + 
(1 — P)62 where /3 is the area fraction of slip length 61. When ^2 ^ ^ L, the effective slip 
length is given by l/6eff = + (1 — /9)/52- These expressions have previously been found 
to hold empirically in molecular dynamics and other numerical simulations of flows over 
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nanostructured superhydrophobic surfaces [l8|. The derivation provided here now provides 
theoretical support for these relationships and elucidates their range of validity. 
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